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ABSTRACT 

We present a map and an angular power spectrum of the anisotropy of the cosmic microwave back- 
ground (CMB) from the first flight of MAXIMA. MAXIMA is a balloon-borne experiment with an array 
of 16 bolometric photometers operated at 100 mK. MAXIMA observed a 124 deg 2 region of the sky 
\ with 10' resolution at frequencies of 150, 240 and 410 GHz. The data were calibrated using in-flight 

> ■ measurements of the CMB dipole anisotropy. A map of the CMB anisotropy was produced from three 

150 and one 240 GHz photometer without need for foreground subtractions. Analysis of this CMB map 
U , yields a power spectrum for the CMB anisotropy over the range 36 < I < 785. The spectrum shows a 

fvi ■ peak with an amplitude of 78 ± 6 f^K at £ ~ 220 and an amplitude varying between ~ 40 (iK and ~ 50 

/iK for 400 < £ < 785. 

m ■ 

> ■ 
m : 

Cn| . 1. INTRODUCTION map the CMB anisotropy over hundreds of square degrees 

t., r , n r,i with an angular resolution of 10'. In this paper we report 

iy-j . Measurements ot the anisotropy ot tne cosmic microwave . f o th MAXIMA 1 fl' ht h'ch t k 1 

background (CMB) can discriminate between cosmological August 2° 1998 6 These include a lOO^s uare°de reeTma 11 
models and determine cosmological parameters with high r U f US ™',n-> •' ese mc u e a square egrees map 

O accuracy (Kamionkowski & Kosowski 1999, and references of the CMB anisotropy and the resulting power spectrum 

> therein). Inflationary dark matter models, for example, over * h « 36 ^ *. f 785 ' ^ch 18 th f lar S est ™ n % e 
,15 ■ j • i c , . i , r reported to date. Despite several common team members, 

. predict a series ol peaks m the angular power spectrum ot ,, , , , • ■ , , r,i , , ,■ , 

5— U f. • , Tn. n . j i. r • the data analysis was independent ot that leading to the rc- 

i ■ the anisotropy. the collected results from many expen- . ^ orted BOOMERANG re Its fdc Bernardi t 1 

O' ments show the existence of a first peak at angular scales ^J^Jf r( ^P or e . * esn s } ?„„!5^ ar n . 1S 6 a " 

L_( ■ !• i ,1 i • i i . ,,. , 2000). A companion paper, Balbi et al. (2000), discusses 

, corresponding to the spherical harmonic multipolc num- . 1 . ... .K, r , „„ , , , 

ber £ ~ 200. These results have been interpreted as evi- thc cosmological significance of the MAXIMA-1 results. 

dence for a flat universe (de Bernardis et al. 2000; Lange et „ . „ 
al. 2000; Dodelson & Knox 2000; Tegmark & Zaldarriaga 

2000; Pierpaoli, Scott, & White 2000). Additional observa- Lee et al. (1999) gives a detailed description of the 

tions probing a broad range of angular scales would greatly MAXIMA system. It is based on a well-baffled, under- 

■ increase confidence in these results and further constrain filled, off-axis Gregorian telescope with a 1.3 m primary 

' cosmological parameters. mirror, mounted on an attitude-controlled balloon-borne 

MAXIMA is a balloon-borne experiment optimized to platform. A well-baffled liquid-Helium-cooled optics box 

is lined with absorbing material (Bock 1994) and con- 
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Table 1 

Central frequencies {u ), frequency bands (f raj „,(/ m «), beam full widths at half maximum (FWHM), detector 

TIME CONSTANTS (t), AND DETECTOR NOISE EQUIVALENT TEMPERATURES (NET) IN CMB THERMODYNAMIC UNITS FOR THE 
FIVE MAXIMA-1 PHOTOMETERS DISCUSSED IN THIS PAPER. THE FWHM, T, AND NET WERE DETERMINED USING FLIGHT 
DATA. THE FREQUENCY BANDS WERE MEASURED IN THE LABORATORY. 



tains two rcimaging mirrors, low-pass filters, field and 
aperture stops, feed horns for the 16 photometers, and 
a focal-plane stimulator. Eight conical single-mode horns 
at 150 GHz and four multi-mode Winston horns each at 
240 GHz and 410 GHz provide 10' beams at all three fre- 
quencies. The frequency bands are defined with absorp- 
tive and metal-mesh filters. Radiation is detected with 
spider-web bolometers (Bock et al. 1996) operated at 0.1 K 
with an adiabatic demagnetization refrigerator (Hagmann 
& Richards 1995). The bolometers are AC-biased to avoid 
low-frequency amplifier noise. Additional channels with a 
constant resistor, a thermometer, and a dark bolometer are 
used to monitor electromagnetic interference, cross-talk, 
and drifts in electronic gain and temperature. The gon- 
dola azimuth is driven by a reaction wheel using informa- 
tion from a two-axis magnetometer and a three-axis rate 
gyroscope. The telescope elevation is set using information 
from an angle encoder. Observations were carried out at 
fixed elevation with the primary mirror scanning ± <~ 2° 
in azimuth at 0.45 Hz and the gondola also scanning in 
azimuth but at a frequency of ~ 0.02 Hz. Both scans were 
triangle functions of time with smoothed turnarounds. 

3. OBSERVATIONS 

The MAXIMA-1 flight was launched from the National 
Scientific Balloon Facility in Palestine, Texas at 1 UT on 
August 2, 1998. Observations of the CMB dipolc for the 
purpose of calibration began at 3.6 UT when the payload 
reached an altitude of 32 km and ended at 4.2 UT after 
~ 100 rotations at 3.3 rpm. The elevation angle was set to 
51°. The payload reached float altitude of 38.4 ± 0.4 km 
at 4.6 UT. 

The 1.6 hour CMB-1 observation began at 4.35 UT with 
a telescope elevation of 46.3°. The gondola was scanned 
±4.1° in azimuth at 16.1 mHz centered at 321.5°. The 
1.4 hour CMB-2 observation began at 6.0 UT with a tele- 
scope elevation of 32.3°. The gondola was scanned ±2.9° 
at 21.3 mHz centered at 323°. Because of sky rotation, the 
combination of these observations covered a nearly square 
region of the sky with an area of 124 square degrees of 
which 45% is cross-linked at an angle of ^22°. 

Observations of Jupiter were carried out from 7.5 to 8.1 
UT to map the telescope beams and provide additional 
calibration information. The elevation was fixed at 44.2° 



while sky rotation and the primary mirror modulation pro- 
vided ~ 200 transits across each beam. 

4. POINTING RECONSTRUCTION AND CALIBRATION 

We identified the stars which moved through the field of 
a CCD camera aligned with the center of the primary mir- 
ror scan by using the balloon location, telescope elevation, 
and the position of Polaris in an offset CCD camera. In- 
terpolations using an angle encoder on the primary mirror, 
rate gyroscopes, and the known star positions permitted 
pointing reconstruction to better than 1' RMS. Less than 
0.1% of the data had pointing uncertainty larger than 2' 
and were not used. 

A full beam calibration of the 150 and 240 GHz pho- 
tometers was obtained from observations of the CMB 
dipole. The data from each rotation were x 2 -htted to a lin- 
ear combination of a dipole model (Lineweaver et al. 1996), 
a galactic-dust emission model (Schlegel, Finkbeiner, & 
Davis 1998), data from one 410 GHz photometer, an off- 
set, and a gradient. The amplitude of each of these com- 
ponents was treated as a free parameter. A monotonic 
change in the detector calibration of less than 9% through- 
out the CMB observations, due to an increasing detec- 
tor temperature, was monitored by illuminating the focal 
plane with the stimulator lamp. Estimated la calibration 
uncertainties were less than 4% for each of the 150 and 240 
GHz photometers. The uncertainties in the dipole cali- 
bration and the time dependent calibration contributed 
about equally to the total error, and systematic sources 
contributed about 25%. 

Beam maps and an independent calibration were ob- 
tained from observations of Jupiter. We estimate a la 
uncertainty of ±0.5' in the size of the beams. The beam 
profiles were integrated and used with the angular diam- 
eter and brightness temperature of Jupiter (Goldin et al. 
1997) and the optical bandpass functions to calibrate all 
16 photometers. For the data reported here, the errors 
in the calibration from Jupiter were between 12 and 14% 
(with about 10, 5, 5, and 2% coming from uncertainties 
in the beam solid angle, frequency bands, Jupiter's flux, 
and measured signal, respectively). The absolute calibra- 
tions from the dipole agreed with those from Jupiter to 
within lcr, with the Jupiter calibration predicting larger 
temperature fluctuations by 11 to 14%. 
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5. MAP AND ANGULAR POWER SPECTRUM 

In this paper we report on the analysis of data from 
the three 150 GHz, one 240 GHz, and one 410 GHz pho- 
tometer described in Table 1. At 150 and 240 GHz these 
photometers are the most sensitive in the MAXIMA array 
and give the highest sensitivity of any CMB instrument 
reported to date. 

The raw data for each photometer consisted of 2.3 mil- 
lion samples of which about 16% were not used. We re- 
moved the stimulator calibration events and other events 
with an amplitude larger than 6a. This procedure broke 
the data into 20 segments that were treated as independent 
observations of the sky. Samples in each of the segments 
which were in excess of 4cr, such as cosmic ray hits and 
short telemetry drop-outs, were removed. For the data of 
one of the photometers, we repeated the data analysis by 
using a threshold of 3a with no significant change in the 
resulting angular power spectrum. We deconvolved the 
transfer functions of the bolometers and readout electron- 
ics and estimated the noise power spectrum from sections 
of the time stream that had no gaps (Stompor et al., in 
preparation). We used the procedure of Ferreira & Jaffe 
(2000) to confirm that the time-domain data are domi- 
nated by noise. We marginalized over frequencies lower 
than 0.1 Hz and higher than 30 Hz, where we did not 
expect appreciable optical signals. 

The calibrated time stream data were combined with 
the pointing solution to produce a maximum likelihood 
pixelized map of temperature anisotropy and a pixel-pixel 
noise correlation matrix for each photometer (Wright 1996; 
Tegmark 1997; Bond ct al. 1999). An area of <~ 20 square 
degrees of the map was not well cross-linked and is not 
included in the present analysis. The data showed a signal 
that was phase-synchronous with the primary mirror scan 
which was due to radio-frequency interference from on- 
board transmitters modulated by the mirror motion. This 
signal was constant within each data segment, varied be- 
tween different photometers, and had an amplitude of 100 
- 300 /xK. We removed it by allocating fictitious map pixels 
to values of the primary mirror angle and determined the 
maximum likelihood map in these pixels simultaneously 
with the temperature anisotropy map (Stompor et al., in 
preparation) . 

We verified that there are no noise correlations between 
maps of different photometers by producing difference 
maps from the data of pairs of photometers. The angular 
power spectra of these maps were consistent with no sig- 
nal. We also showed that histograms of the temperatures 
in the difference maps were consistent with the distribu- 
tions expected for no noise correlations at a Kolmogorov- 
Smirnov significance level larger than 10%. A combined 
temperature anisotropy map was then produced by adding 
individual maps with a weight inversely proportional to 
their noise correlation matrices. A Wiener filtered version 
of this map is shown in Figure 1. We assign a calibra- 
tion uncertainty of 4% to the magnitude of temperature 
fluctuations in the combined map. 

We calculated the angular power spectrum Ce of the 
combined map using the MADCAP (Borrill 1999) imple- 
mentation of the maximization of the likelihood function 
following Bond, Jaffe & Knox (1998). This implemen- 
tation assumes that the beam shape has axial symmetry. 



We produced an effective beam for the analysis of the com- 
bined map by noise- weight averaging the individual beams. 
We followed the procedure of Wu et al. (2000) to find a 
symmetric approximation for the effective beam and in- 
cluded the small smoothing provided by the pixelization. 
We tested this procedure for the MAXIMA- 1 beams and 
data and found no systematic bias of the C( estimates 
(Wu ct al. 2000) . We calculated the power spectrum 
of the temperature fluctuations using 15 bins in I over 
the range 3 < £ < 1500 assuming a constant £(£ + l)Ci 
band power in each bin, and marginalizing over the CMB 
monopole and dipole. We further marginalized over the 
bins at £ < 35 and £ > 785 and diagonalized the £-bin cor- 
relation matrix using a variant of techniques discussed in 
Bond, Jaffe & Knox (1998). The correlations between the 
dominant bin and adjacent bins were typically less than 
10%. Table 2 lists the dominant bins, the Ce estimates, 
and the AT = \J £(£ + l)Ci/2ir estimates for the corre- 
sponding uncorrelated linear combinations of bins. We 
quote la errors on the Ce estimates assuming 68% con- 
fidence intervals using the offset log-normal distribution 
model of Bond, Jaffe & Knox (2000). These errors do not 
include two additional independent sources of uncertainty. 
Expressed as a percentage of £(£ + l)Ce/2n, the la calibra- 
tion error is a constant 8% and the ^-dependent error due 
to the beam-diameter uncertainty, which has been shown 
to be fully correlated between £ bins, is given in Table 2. 
Information on the shape of the bin-power likelihood func- 
tions and window functions will be made available on the 
MAXIMA web site (http://cfpa.berkeley.edu/maxima). 

The top panel of Figure 2 shows the maximum likelihood 
power spectrum, an inflationary adiabatic model that best 
fits the MAXIMA- 1 and COBE/DMR power spectra, and 
a ACDM model. The x 2 for the best fit and for the ACDM 
models are 41 and 50, respectively, using all 38 data points. 
If we only use the 10 data points of MAXIMA-1 we ob- 
tain x 2 = 7 and 9 for the best fit and ACDM, respectively 
(Balbi et al. 2000). 

6. FOREGROUNDS 

Foreground sources include emission from the earth, 
the atmosphere, interstellar dust, free- free radiation, syn- 
chrotron radiation, and point sources, and scattering due 
to the Sunyaev-Zeldovich effect. The two CMB scans, per- 
formed with a separation of 1.5 hours at different telescope 
elevations, show consistent structure. This temporal sta- 
bility is inconsistent with an atmospheric or ground-based 
origin for the signal. We extrapolated the 100 yum, 10' 
resolution Schlegel, Finkbcincr, & Davis (1998) map of 
our observed region to lower frequencies using the two- 
component dust model favored by Finkbeiner, Davis & 
Schlegel (1999). The predicted RMS dust temperature is 
only 2.3 and 9.3 fj,K, for the 150 and 240 GHz bands, 
respectively. We found a statistically significant corre- 
lation between the dust model and our maps. The ra- 
tio of the detected to predicted RMS dust signal was 
statistically consistent with unity. When we subtracted 
the correlated dust signal from the maps the change in 
the measured angular power spectrum was negligible. A 
catalog search (Sokasian, Gawiser, & Smoot 2000; Ga- 
wiser & Smoot 1997) yielded no detectable radio or infra- 
red sources in any of the frequency bands. Estimates 
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Table 2 

Uncorrelated Power Spectrum from the MAXIMA-1 map. Errors are the 68% integrated pro bability of t he 

OFFSET LOG-NORMAL LIKELIHOOD FUNCTIONS WITH A CONSTANT PRIOR IN EITHER l[l + l)C e /2lT OR AT = \J £(£ + l)Ci/2lV FOR 
THEIR RESPECTIVE COLUMNS. HERE [£ m in,£max] GIVES THE RANGES OF THE DOMINANT BINS. THE CORRELATED BEAM ERRORS 

ARE DESCRIBED IN THE TEXT. 



of brcmsstrahlung and synchrotron radiation (Bouchet & 
Gispert 1999) yielded contributions of less than 1 /xK at 
150 and 240 GHz and no subtractions were made. 

Integrating the measured power spectrum of one pho- 
tometer at 150 GHz and the one at 240 GHz we find 
for a thermodynamic temperature fluctuations ratio of 
0.91 ± 0.18 compared to unity for the CMB, 0.06 for emis- 
sion from dust, 2.1 for synchrotron, and 3.7 for free- free 
radiation, assuming the Tegmark et al. (2000) "middle" 
foregrounds model. For the ratio between 150 and 410 
GHz, we found a lower limit of 1.1 • 10~ 2 compared to 
3.2 • 10~ 4 for dust and 5.7 • 10" 5 for atmosphere. The ob- 
served anisotropy is not consistent with the thermal SZ 
effect, which would produce anti-correlated structure in 
the 150 and 240 GHz bands. 

7. TESTS FOR SYSTEMATIC ERRORS 

Because of computational limitations, tests for system- 
atic errors were done with maps which had square pixels 
of 8' and 10' on a side. The power spectra calculated from 
the 8'-pixel and 5'-pixel combined maps were statistically 
consistent. The following combinations of the data were 
analyzed and produced a power spectrum consistent with 
no signal: (1) a dark bolometer, (2) the data from the 410 
GHz photometer, (3) the difference between the overlap- 
ping part of the combined map from the CMB-1 and -2 
scans, (4) the differences between the maps produced by 
different photometers. We also weight-averaged the maps 
of the second and third photometers in Table 1 and the 
first and fourth. The maximum likelihood estimate of the 
power spectrum of the difference between these indepen- 
dent maps is consistent with no signal as shown by the 
open circles in the top panel of Figure 2. 

We compared the estimate of the angular power spec- 
trum to that obtained using: (1) only the sections of the 
map where the CMB-1 and -2 scans overlap, (2) a map of 
each of the CMB scans alone. We made maps and calcu- 
lated Ci estimates from the data of each photometer alone 
and using: (1) only a sub-section of the time stream data, 
(2) a high-pass filtered version of the time stream where 



the high-pass was a time domain box-car with a width of 
10 sec, (3) various combinations of frequency marginal- 
izations between 30 and 70 Hz, and 0.05 and 0.3 Hz, re- 
spectively. In all these cases the computed power spectra 
agreed among themselves and with the power spectrum 
presented in this paper. The Kolmogorov-Smirnov test 
(see Section 5), as applied to difference maps and to a 
map of the dark bolometer, confirms that we correctly 
estimated the noise in the experiment and that the pixel- 
domain noise is Gaussian. 

We used simulations to test the algorithm used to sub- 
tract the signal that was phase-synchronous with the pri- 
mary mirror modulation. We found that the power spec- 
trum estimate was not biased with phase-synchronous sig- 
nals larger than those observed in the data. If we make 
the maps without removing the phase-synchronous signal 
the power spectrum estimate changes only at £ < 200. 

The computer programs used to generate the maps and 
power spectra were tested extensively using simulations of 
the time domain data and noise. Maps were produced by 
two independent computer codes and the power spectra 
calculated from these maps were consistent. We used one 
of the map-making codes to make maps of Jupiter and 
found them consistent with those obtained with a simple 
data-binning technique. 

8. DISCUSSION 

We have observed temperature anisotropy on the sky at 
150 and 240 GHz that is consistent with fluctuations in 
the cosmic microwave background radiation, and inconsis- 
tent with any known foreground. The observations were 
carried out with photometers that give the highest CMB 
sensitivity reported to date. Our measurements cover a 
range of angular scales corresponding to the multipolc 
range 36 < £ < 785, which is the largest yet reported by 
a single experiment. The measured angular power spec- 
trum shows a clear peak at £ c± 220, and an amplitude 
varying between ~ 40 fiK and ~ 50 fiK at 400 < £ < 785. 
The power spectrum is well fit by an inflationary adia- 
batic model over the entire range of £. The best-fit model 
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has a total energy density close to unity and a non-zero 
cosmological constant. The MAXIMA- 1 power spectrum 
appears consistent with that of the BOOMERANG exper- 
iment (de Bernardis et al. 2000) once the power spectra 
of the two experiments are scaled by factors equal to their 
respective la calibration uncertainties, see Figure 2. A de- 
tailed analysis of the combined data sets which includes a 
determination of the calibration factors that bring the ex- 
periments to agreement is presented in Jaffe et al. (2000). 
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Fig. 1. — A Wiener filtered map of the CMB anisotropy with 10' resolution from MAXIMA-1. The map is made using data from three 150 
GHz and one 240 GHz photometer and contains 15,000 5' x5' pixels. We used the angular power spectrum shown in the top panel of Figure 2 
as the prior for the Wiener filter. Pixel boundaries have been smoothed using interpolation. 
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Fig. 2. — Top panel: Angular power spectrum of the CMB anisotropy from the MAXIMA-1 map shown, in Figure 1 (filled circles). The error 
bars are 68% confidence intervals calculated using the offset log- normal likelihood functions of Bond, Jaffe & Knox (2000). The solid (dashed) 
curve is the best fit (ACDM) inflationary adiabatic cosmology to the MAXIMA-f and COBE/DMR data. The models have (f2f,, f2 ct j m , SIa, 
n, h) = (O.f, 0.6, 0.3, 1.08, 0.53), and (0.05, 0.35, 0.6,1 ,0.65), respectively (Balbi et al. 2000). The open circles are the estimated power 
spectrum of the difference between two independent maps, each produced by a weighted-average of the maps from a pair of photometers. 
Bottom panel: A comparison of the MAXIMA power spectrum with that of the recently reported BOOMERANG experiment (de Bernardis 
et al. 2000). Consistency between the power spectra has been achieved by scaling the MAXIMA-1 (BOOMERANG) power spectrum down 
(up) by a factor equal to its 8% (20%) lcr calibration uncertainty. These data show a suggestion of a peak at i ~ 525. 



